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Abstract 


The device considered in this st’ady is a semiconductor waveguide 
isolator consisting of a hollow colurrin of a semiconductor mounted coaxially 
in a circular v?aveguide in a iongitudinal dc magnetic field. An elementary 
and physical analysis based on the excitation of plans waves in the guide 
and a more rigorous mode raatching analysis are presented. These theoretical 
predictions are compared with experimental results for an InSb isolator at 


94GHz and 75 °K. 
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I, INTRODUCTION 

The purpose of this paper is to present a theoretical study of a 

semiconductor waveguide isolator of new geometry and the experimental 

results for an InSb isolator at 94 GHz and 75°K. 

The propagation characteristics of em waves in a circular waveguide 

partially filled with a plasma have been investigated extensively examining 

[ 1 ] 

the effect of a plasma upon waveguide modes. For example^ R. N. Ca.rlx.le ■■ 
measured the propagation constants of the modes in a cylindrical vrave- 

guide partially filled with a gaseous plasma in a longitudinal dc magnetic 
field and observed that right and left circularly polarized (cp) waves 
travel with different propagation constants. The theoretical and experr- 
mental study of a circular 35 GHz waveguide containing a coaxial TnSb rod 
was done extensively by B. R. McLeod and W. G. May'^ ' ' 

The configuration to be studied here, shown in Fig. 1, has a hollow 
plasma column in a circular waveguide. This geometry is more practical 
to construct than the rod configuration for millimeter v^avaguides. The 
theoretical analysis will be in two parts: First, a heuristic argiiment 

based on the excitation of plane waves in the guide is presented. Second, 
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a more rigorous mode matching analysis (henceforth abbreviated MMA) is 
given . 


II. FORMULATION OF THE PROBLEM 

With a longitudinal dc magnetic field on the plasma column, the total 
dielectric constant of the material seen by the em waves is due to both 
lattice and free carriers and depends on the direction of the circular 
polarization (cp) with respect to the longitudinal magnetic field. The 
total relative dielectric constant has the form^^^ : 


= Kj^ ± jK^ = + jk" , 
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Here K is the relative static lattice dielectric constant of the magneto- 
L 

plasma, and t is the carrier relaxation time. The subscript ± refers to 
the left and right cp components in the direction of the dc magnetic field. 
For suitable choice of parameters, could be relatively large and K| very 
low or even negative, and K" small. This possibility is put to use in an 
isolator as described belov; by exciting a right or left cp wave in the 
circular waveguide containing an annular plasma column in a longitudinal 
magnetic field, as shown in Fig. 1. 
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II-l. APPROXIMATE THEORY FOR ATTENUATION COEFFICIENT 

In this section we develop a simple physical picture from which the 
performance of the structure of Fig. 1 can be estimated. Using a heuristic 
argument, we will find the ratio of the electric field strength in the 
annular plasma column to the field in the center section of air-filled 
region. It is assumed that the fields both in the center section of the 
air-filled region and in the annular plasma column are plane waves of 
approximately the mode pattern. The amplitudes of fields in the two 

regions can be, however, drastically different and depend on the dielectric 
constant of each region. We further assume that a simple static analysis 
will suffice to determine these relative amplitudes. Of course, the above 
assumptions are not truly valid and will not be used i.n the more rigorous 
MMA in the next section. 

The sim.ple. static analysis gives for the ratio of the field concen- 
trations ; 
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1 + 
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where K is the ratio of the relative dielectric constants of the outer and 
inner media, and E^ and are the "uniform" fields in the inner and outer 
regions. Since the inner region is empty, K will be interpreted as K^, 

Eq. 1, for cp waves. 

If the material parameters and magnetic field are such that 
1 /t and if the dc magnetic field is adjusted so that as in (1) is close 
to zero (at least K|<<K^), it appears reasonable from (6) that the left cp 
waves are strong3.y excited in the outer region, the annular plasma column. 
Thus f the wave propagates almost exclusively through the p2.asma, and be-- 
cause of tlie lossy part of the dielectric constant of the plasma , much of 
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the left cp wave will be absorbed by the plasma and dissipated as heat. On 
the other hand, the same era wave traveling in the opposite direction wrth 
respect to the longitudinal dc magnetic field is right cp. It is predicted 
from (6) that this wave is largely excluded from the annular plasma column 
because K' will be large for InSb, 32) , and with the wave will 

be only slightly absorbed by the plasma. 

Let us examine the attenuation coefficient using the approxiiaate model 
discussed above. Consider an infinitely long annular" plasma column placed 
in a uniform era field with the propagation direction parallel to the axis 
of the column. A.ssuming that the transverse fields in the annular plasm.a 
column are aporoximately uniform for both polarisations of the era vraves, 
the power dissipation per unit length in the plasma annulus is 
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where E is the field inside the annular plasma column, A is its cross- 
o 

sectional area and is the effective conductivity of the plasmaj, ue^K" . 
From the foregoing discussion, the field inside the plasma column is also 
assumed to be a uniform TE mode, and the power flow through the plasma 
column, P^, is then approximately given by 
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where is the impedance of the TE mode in the column, - 377 ohms, and 
X is the guide cutoff wavelength. Similarly the power flow through the 


center air-filled region, P^, is given crudely by 

1 l^ii' 
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where E. is the field in this region and A^ is the area of its cross section 
The approximate attenuation coefficient, a^, due to absorption in the plasma 
found by dividing the dissipated power (7) by the total power flow 
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and is found to be: 
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This expression will be evaluated and compared to the more accurate 
expression derived from the MMA in the next section. 


I I -2. MODS MATCHING ANALYSIS FOR THE BOUNDARY VALUE PROBLEM 

The general calculations of the reflection and transmission coeffi- 
cients of a finite-length section of the annular magnetoplasma column as 
in Fig. 1 consist of the following two steps: 

.1.) Determi.nation of the propagation constants of the 
gyromagnetic modes in the section of the annular 
magnetoplasma column, 0<z^S,, 

2) Use of these propagation constants to match the 
boundary conditions at the ends of plasma column, 
z=0,£. 

The radial bcondary conditions are that the tangential components of E are 
zero at r=b and that the tangential components of E and H are continuous 
at r=a. The characteristic equation turns out to be a 6x6 determinant (see 
discussion in Appendix A) . The problem of matching fields at the ends of 
the plasma column is quite complicated since a single incident mode can 
excite a number of gyromagnetic modes in the plasma section and these, in 
turn, excite a number of TE and TM. modes in the empty waveguide, and these 
higher modes must be included in applying the boundary conditions. Further- 
more, the gyromagnetic modes do not form orthogonal sets due to the existence 
of .loss in the plasma [6], and we use the point matching method to compute 



approximately the reflection and transmission coefficients at the ends of 
the plasma column, by matching boundary conditions at only a finite number 
of points over the cross section. This formulation becomes increasingly 
exact by increasing the number of modes considered and the number of 
matching points chosen. The overall attenuation due to the finite section 
of the annular plasma column is then calculated from the transmission 
coefficient of the dominant mode at the end of the plasma column. 

The details of the calculations are given in Appendix B, and some 
numberical results v/ill be presented in §IV below. 


III. experimental procedure 

The material chosen for the solid state plasma "was n-type InSb. This 
material had n = 2 . exlO^'^^/cm^ and y = 5 . 2xl0^cm^/V-sec at the operating 
temperature (VS'^K) , and we assumed that - 16.0 and m* = 0.014m^. 

The isolators tested consisted of two transitions and the test section 
shown in Fig. 2. The transitions convert the TE^^ mode propagating in the 

rectangular waveguide into the cp the circular waveguide. The 

test section contains the annular column of n-type InSb, which was cut by 

using an ultra-sonic drill and polished chemically in a bromine and methanol 
etch. The small opening between the column and the circular waveguide wall 
was filled with highly conducting silver paint, and styrofoam filled the 
center section. The devices were placed in a longitudinal dc magnetic 
field using a 6-inch Varian magnet with a 2-inch gap. A styrofoam Dewar 
surrounding the semiconductor was filled with the liquid nitrogen which was 
not allov.'ed to enter the test section. The microwave system was conventional 
with a precision attenuator. A block diagram of the experimental set-up is 
shown in Fig. 3. 
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IV-2. COMPARISON WITH MODE MATCHING ANALYSIS 

The first step of MMA was to determine the propagation constants of 
the gyromagnetic inodes. Since the characteristic equation is quite com- 
plicated, an iterative root-finding technique was employed to find the 
propagation constants. The real and imaginary parts of the propagation 
constants for both right and left polarizations were obtained by computer 
to three place accuracy. 

The resulting propagation constants were utilized to match the 
boundary conditions at the ends of the plasma column. This second step 
determines the amplitudes of the reflection and transmission coefficients 
of the gyromagnetic modes and of the empty waveguide modes. In the present 
case, four gyromagnetic modes in the section of the InSb column and tvio TE 
and two TM modes in the empty waveguide were taken into account. Point 
matching was employed at r=0.3 and 0.7mrri to solve the inhomogeneous linear 
equations of order 16. The overall attenuation due to the various lengths 
of the InSb columns v/as evaluated from the transmission coefficients of TE, ^ 

j. 

circular mode at the ends of the column and is shown in Figs. 4 and 5 by 
the solid lines . The general agreement between the experimental and the 
theoretical results is quite good. The boundary conditions are imposed at 
only four points at the ends of the plasma column by taking into account 

«r 

only four gyromagnetic modes in the raagnetoplasma section and four higher 
excited modes in the empty waveguide section. It is, however, expected 
that, although the inclusion of higher excited modes as well as the dominant 
modes is essential in a rigorous treatment of the boundary conditions at 
the ends of the plasma column, the convergence is fairly rapid [7] . Thus 
small discrepancy between theoretical and experimental results appear to 
be mainly due to lack of perfect cylindrical symmetry in the experimental 
setup. 
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APPENDIX A 

CHAPACTERISTIC EQUATION FOR CIRCULAR 
WAVEGUIDE CONTAINING ANNULAR PLASMA COLUMN 


This appendix outlines the boundary value problem that results from 
the configuration shov;n in Fig. 1. 

The appropriate solutions of the v/ave equations for the longitudinal 
field components for r<a are 

E° -- UJ (qr) , (11) 

z m 

and 


11° = VJ. , (12) 

where U and V are arbitrary, complex constants and F the complex propaga- 
tion constant. Corresponding transverse E field components for r<a are 

' [-"I , (13) 


= T-u-S 


lf^„(qr) . .;(qr)] 

q r U- 


with similar expressions for H. 

In the region of the annular plasma column a<r<b. Maxwell's equations 
yield the coupled wave equations: 
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where a, b, c, d can be defined in terms of F and the plasma dielectric 
tensor components. Assuming solutions of the form exp(-jp*r,p) , these can 
be solved giving the dispersion equation: 
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- (a+c)p^ + (ac“bd) = 0 , (17) 

The longitudinal E field inside the annular plasma (a<r<b) is 

^2 " DN^(p2r) ] , (IS) 

with a similar expression for H , and A, B, C and D are arbitrary^ com- 
plex constants with similar but more complicated expressions for the 
other components of E and H, 

Now the boundary conditions at the plasma surface (r=a) shall be 
taken as the continuity of the tangential E and H fields; (4 equations) . 
Also at thie boundary of the perfect conductor waveguide wall (r-b) , tangen- 
tial E and H fields are zero; (2 more equations) . In order to have non- 
trivial solutions of A, B, C, D, U and V, the 6x6 determinant of these 
coefficients for the 6 equations must be zero. The details will not be 

published here due to space limitations but v^ill be supplied upon request 

[ 8 ] 

to the authors 
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APPENDIX B 

I4ATCHIWG FIELDS AT THE ENDS OF THE 
MAGNETOPLASMA COLUMN 


For a plasma section of finite length, one is confronted with the 
problem of matching fields at the ends of the plasma column. Let us con- 
sider the configuration shown in Fig. 1 when the incident wave is TEj^^. 


One coiTiT:)onent of the H field in region 1 is 

jmcf)-r®z g 
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jm(fi+r^z 


(19) 


where R , is a reflection coerficient, with similar expressions for the 
ni 

„E 

other components of H and for E. The mode propagation constant, is 


determined from the characteristic equation of the TE wave, i.e., 


where 
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In the plasma a<r<b, one component of E field is 
n L 

-1 z j- 
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and for o<r<a 
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" . . for the other field components of H and . 

Lth corresponding expressi efficients for the n'*^^ TE (or j-M) 

The reflection and transmission co- ^^termined by the 
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FIGURE CAPTIONS 


Figure 1 . Circular Waveguide Containing an Annular InSb Column 
Figure 2. Waveguide Isolator Containing Annular Plasma Column 
Figure 3 . Block Diagram of Experiment Set-up 

Figure 4 . Theoretical and Experimental Results for Hollov/ InSb 
Waveguide Isolator 

[Geometry of 94 GHz InSb vJaveguide Isolator is shown 
in Fig. 2 with •t=0. 762mm.] 


Figure 5. Theoretical and Experimental Results for Hollow InSb 
Waveguide Isolator 

[Geometry of 94 GHz InSb Waveguide Isolator is shown 
in Fig. 2 v/ith -1-1.5241™.] 
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